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Abstract 

Based on a perusal of available experimental data, we propose a possible 
explanation for the heavy-fermion Fermi liquid state coexisting with quasi- 
one dimensional antiferromagnetism in the low-carrier pnictide Yb^As^. Both 
of these low temperature features are shown to originate from one single fact 
- the charge ordered state forming at high-T is imperfect, due to quantum 
fluctuations in the CO state. The quasi-lD magnetism is understood in terms 
of the physics of the squeezed Heisenberg chain, while the heavy fermion-like 
behavior arises from the compensation of the fraction of the Yb moments 
(which are disordered) by the 2p band of As. The proposed picture shows 
how the seemingly strange features observed experimentally can be reconciled 

as a result of this single hypothesis. 
PACS numbers: 71.28+d,71.30+h,72.10-d 
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1. INTRODUCTION 

Rare-earth pnictides like Yb 4 As^ have been increasingly studied in the past few years, as 
they are thought of as being " low-carrier density Kondo systems" . Transport measurements 
show similarities with typical heavy fermion Fermi liquid metals [1], while the thermody- 
namics and the magnetic response [2] is characteristic of one-dimensional spin systems. A 
naive estimation of the carrier number from the Hall data leads one to conclude that the 
carrier density is very low [3]. How such a supposedly small carrier density can screen the 
Yb local moments to give a heavy FL at low-T is the issue central to the understanding of 
this material. Purely on theoretical grounds, one might want to ask: 

(1) Under what conditions can a heavy-fermion like FL metal coexist with quasi-lD spin 
fluctuations? How do these degrees of freedom affect each other at low-T? 

(2) A more general, but related question: Are there hitherto unknown routes to heavy 
fermion behavior, other than the well-studied Kondo effect? 

In this letter, we want to propose a route which ties together the various seemingly 
irrevocable facts that are known about this material. We start with a brief recapitulation 
of known experimental and theoretical facts, and then go on to show how both the heavy- 
fermion like transport and the lower- dimensional magnetism arise naturally from imperfect 
nature of the charge ordering that occurs at higher T. Our analysis also provides a natural 
explanation of the fact that the closely related materials Yb 4 X 3 ,X = P, Sb, behave very 
differently from Yb^As^. The phosphide is an insulator all the way down to the lowest T, 
with a more perfect charge order, while the antimonide is a mixed-valence metal (no charge 
ordering) down to the lowest T. 

Yb 4 As3 is a semiconductor at room temperature, with the rare-earth /-ion of Yb in a 
mixed valence state {Yb 1A+ — Yb 13+ ). At about 293 K, a collective JT distortion causes 
three of the four identical diagonals along the cubic cell to elongate, leading to a contraction 
of the fourth diagonal, into which the Yb 13+ ions gather. Since Yb 13+ carries a spin 1/2, 
one has an (almost half-filled) chain with strong electronic correlations. These are precisely 
the degrees of freedom giving rise to the ID spin fluctuations observed in inelastic neutron 
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scattering, and to the spin susceptibility and the 7 coefficient of the specific heat, as shown 
on [2], and discussed on experimental [2] and theoretical grounds [3]. 

Since the spin degrees of freedom are those of the ID Hubbard model, the spin excitations 
are fermionic spinons rather than conventional bosonic spin waves. But this is not all, as 
bandstructure calculations reveal that the carriers seem to come predominantly from the As 
2p-band. A consistent picture of the low-T transport must necessarily involve the scattering 
of these carriers on the spin fluctuations of the Hubbard chain. If the chains were half- 
filled, leaving the As p carriers out of the reasoning would lead one to conclude that the 
material must be a Mott-Hubbard insulator. On the other hand, if one does include the As 
carriers, a heavy fermion-like resistivity can only result from the scattering of these light 
carriers off almost localized Yb 4f-spins. The issue then is whether this scattering (details 
have not been worked out) can give rise to a quadratic (in T) scattering rate. The large A 
value suggests strong scattering, but in view of the fact that one is dealing with spinon-pair 
fluctuations that scatter the light carriers, a concrete calculation is needed to answer the 
question. We argue, however, (see below) that this mechanism is not responsible for the 
observed resistivity, which, we show, has a more conventional interpretation. 

This is intimately related to the second question posed above: the scattering giving rise 
to a HF-FL metal at low-T must necessarily involve quenching the local moment degree of 
freedom of the Hubbard chain. So one might want to address the question of a Kondo-like 
effect originating in a way different from the usual accepted one. The basic question should 
then be rephrased as: 

Can light electron-like carriers scattering on spinon-pair fluctuations yield heavy fermion 
like behavior? 
or, 

Is the above view incorrect, and could a conventional Kondo effect occur as a result of 
scattering of the p-holes off the spin fluctuations of the fraction of the localized spins in the 
other three chains? 

If the second point of view indeed holds, a consistent explanation for the above anomalous 
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features can indeed be given without having to take recourse to exotic scenarios. 

To be able to decide the issue, one has to appeal to experiment. The experimental results 
need to be correlated and a common thread connecting the observed phenomena needs to 
be searched. We summarize the experimental results below. 

(1) The high-T phase is a mixed- valent metal. This undergoes a first-order charge- 
ordering transition at T c = 293K via a collective band-JT effect. The effect is to shorten 
one of the four equivalent chains (at high T) and to elongate the other three. The /-holes 
carrying a spin S = 1/2 sit on the shorter chains, and because the Hubbard Uff » t/, the 
effective / — / hopping strength, one has to deal with a strongly correlated ID Hubbard like 
chain. But this is not the whole story, because bandstructure calculations have shown that 
the carriers responsible for conduction are in fact derived from the As 2p-band (which is a 
wide band, giving rise to small effective- mass carriers). The low-T heavy fermion behavior 
has to do with the scattering of the light p-holes off the almost localized /-moments of the 
Yb. The question for theory is: How can this give rise to a resistivity, p(T) ~ AT 2 with 
large Al 

(2) The specific heat shows a linear T-dependence, and one might naively associate this 
as one of the signatures of usual heavy- fermion behavior especially since the ratio A/7 2 is 
also characteristic of heavy fermion metals. What gives the game away is that a related 
material, F6 4 P 3 , an insulator, also has the linear term, with a 7-coefficient larger than that 
for Yb 4 As 3 \ Thus, the linear specific heat can only result from the spin degrees of freedom, 
spinons, of the ID Hubbard model. That this is indeed the case is shown conclusively by 
inelastic neutron scattering data, which reveals a spin dispersion which is almost that of 
an ideal ID S — 1/2 AFM chain. An external magnetic field has been shown to result in 
significant changes in the 7 coefficient, as well as in the spin susceptibility, while leaving 
the resistivity almost unchanged! Exactly the opposite happens with the effect of external 
pressure, as we discuss below. 

First, we summarize the effect of an external magnetic field. A uniform longitudinal field 
has been shown to generate a staggered transverse field component along the x-direction, 
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resulting in the effective model given by 

H ch = 2lJ2 Si.S, + i - h £(-l)<5? (1) 

i i i 

Calculation of the specific heat, C(T), shows the opening up of a gap at low T in good 
agreement with experimental observations [4]. The static magnetic susceptibility depends 
on the orientation of the external field relative to the spin chain. What is so surprising 
is that this change in the spin excitation spectrum has very little influence on the low-T 
resistivity. This is hard to understand if the ID spin excitations are assumed to drive the 
observed Fermi liquid resistivity. 

Application of hydrostatic pressure has exactly the reverse effect: the resistivity changes 
drastically as pressure is increased, the main features being the suppression of the high-T 
"charge ordering" transition, a gradual but strong reduction of the broad maximum in p(T), 
and a drastic decrease in the coefficient A. The suppression of the high-T CO state can be 
understood in terms of the increased hybridization between the As p and the Yb 4/ orbitals. 
The specific heat and spin susceptibility, however, are completely unaffected by pressure. 

Replacement of As by an ion like P should stabilize the charge ordered (CO) state. This 
is indeed observed [5], but with a larger 7 coefficient of the specific heat than in Yb 4 As 3 . The 
resistivity in this case is insulating all the way down to low T, confirming that the HF-FL 
behavior in Yb^As?, is not related to scattering off the ID spinons, but is connected with the 
imperfect CO state in Yb 4 As 3 . Experimental estimates [6] point to the disordering of the 
CO state with a small fraction of the Yb 3+ on the three longer chains. This observation is 
further supported by LDA+U calculations [6], which show that the hopping integral between 
the short chain and the longer chains is greater than the intrachain hopping integral in the 
short chain. These S = 1/2 spins are randomly distributed in the three-dimensional host, 
and scatter the light As holes. 

In this paper, we propose that both, (i) the one-dimensional magnetism , and (ii) the 
low-T heavy-fermion FL metallic phase observed in Yb^As^, result from one single phe- 
nomenon: the incomplete nature of the the CO state, which sets in as a first-order transition 
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at T ~ 293K [1]. This high-T CO state has been sought to be explained by a collective 
band Jahn- Teller distortion [7,8], whereby the homogeneous mixed- valent state with non- 
integral Yb valence state of 2.25 is unstable to a state where the Yb 3+ (S = 1/2) congregate 
predominantly on one chain. This chain is shortened, while the other three are elongated. 
However, this CO is not perfect and estimates as well as an explicit calculation of a Hub- 
bard model with the Jahn- Teller distortions in the geometry of four interpenetrating chains 
relevant to Yb^As^ [6] show that a small fraction of the Yb 3+ sites are on the other chains, 
an observation that is reconcilable with the fact that the hopping integral between the short 
and the long chains is greater than that along the short chain. This observation is crucial 
for the picture we propose below. 

The ID magnetism observed in this material is then understood in terms of the physics 
of the ID Hubbard (or squeezed Heisenberg model [9]) chain. As is known [10], the spectrum 
of the ID Hubbard model is characterized by spin-charge separation. The spin is carried by 
neutral fermionic excitations, i.e, by the spinons, while the charge degrees of freedom are 
described by spinless, hard core bosons, the holons. Because of spin-charge separation, one 
can continue using the Heisenberg model (with an exchange coupling that depends on the 
deviation from half-filling) to describe the spin dynamics [11]. 

We have computed the uniform static spin susceptibility and the specific heat from the 
thermodynamic Bethe ansatz (TBA) equations [12] for a modified Heisenberg chain with an 
exchange strength (dependent on the hole concentration) given by 



We point out that such an approach has also been used by Zvyagin [17] in the context 
of Yb^As^. However, while we agree with [17] on the issue of the quasi-one dimensional 
magnetism in Yb^As^, we differ on the origin of the FL-like resistivity (see below). The 
Hamiltonian is, 



I e ff(n) 



n[l — sin(27m)/27m] 



(2) 




(3) 
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It is well known that the ID Heisenberg spin model is completely integrable, and that 
the TBA gives an exact expression for the free energy of the system. Using the spinon 
formulation of the TBA (using the quantum transfer matrix) [12] the problem reduces to 
the solution of two coupled nonlinear integral equations. 

lno(x) = - * Ief f +(k* ln9t)(ar) — (k * lnW)(x - in + ie), (4) 
T cosh x 

lna(x) = - ***** +(k* \nW)(x) - (k * In 91) (a; + in - ie), 

J- COSH Ou 

2l(x) = 1 + a(x), %{x) = 1 + a(x). 

1 r°° e -f fe 

k(x) := — / — — cos(kx)dk, (5) 

v ; 2tt J 2 cosh f k V ; V ; 

where the asterik denotes convolution (/ * <7)(x) := f^oJi 00 ~ y)d(y) ( i x - The free energy is 
given by /(0) with 

= _ r /M« d!/i 

7 cosh(x-y) 

ignoring the contribution of the groundstate energy. Solving these equations by the fast- 
fourier transform method, we have computed the spin susceptibility and the specific heat as 
a function of T for the case applicable to Yb^As?,, i.e, with a small fraction of spins in the 
other three chains (n = 0.9). Calculation of the staggered spin susceptibility has already 
been carried out in [4], and is consistent with a gapped spin- wave spectrum resulting from 
a field-induced anisotropy. 

The calculated curves are in fair agreement with the experimental results for this ma- 
terial. In this context, it is worthy of mention that the effect of static disorder induced by 
the vacancies (absent spins) on the shorter chain may become important at very low T [13]. 
This effect is absent in the above calculation. 

To account for disorder, we use the results in [16] where an arbitrary distribution of 
impurities (which do not destroy the integrability) parametrized by the real numbers 9j at 
the site j on the spin chain, are added to a Heisenberg spin chain. This approach has been 
applied by Zvyagin [17] for the same purpose in the case of Yb^As^ recently. The total free 



energy of the spin chain with impurities is F = L" 1 J2j /(f @j) where the sum is taken over 
all the sites (for sites without impurities we take 9j = 0), and f(0j) is given in (6). Now 
an appropiate random distribution of values 9j, P(9j) oc exp(— 7r(A — can be taken 

resulting in a power law divergence of the susceptibility at very low temperatures and the 
linear coefficient of the specific heat 7 

( X ) oc (7) ~ T A -\ 

taken A ~ 0.26 — 0.42 coincide with data on real disordered quasi ID spin 1 AF systems. This 
is quite consistent with the observed power-law divergence in (%). It would be interesting 
to see whether (7) also shows a similar divergence at very low T. 

The specific heat shows the linear low-T contribution. A comparison of the curves for 
two values of n shows that the 7 coefficient of the linear term increases as n increases towards 
unity, providing a simple explanation for the larger 7 value for compared to Yb^As^. 

The specific heat in an external magnetic field has been studied [4] by mapping the spin 
model in a staggered transverse field (induced by the uniform longitudinal field) onto a sine- 
Gordon model, and accounts naturally for the field-induced opening up of the gap in this 
material. Thus, the magnetic fluctuations are well understood in terms of the excitation 
spectrum of the ID Heisenberg model, as was inferred previously [1-4]. 

To understand the origin of the Fermi liquid characteristics in transport , we observe 
that the small fraction of the localized spins in the three elongated chains can be modeled 
by an Anderson-like model with large U (since the Yb holes in the other three chains do not 
order). The Hamiltonian is then given by, 

H = Y1 C ka C ka + Yl tk ^t C ka + h.c) + ^ Un itf n ifl + J^^icT + E f n iffJ ) (7) 
ka ika i ia 

where the f ia , j\ a represent the Yb 4f fermions localized on the longer chains. In this case, we 
use the approach developed in [14] to compute the local spectral density and the resistivity 
in d — oo. In this case, our effective Anderson model yields a correlated Fermi liquid 
metallic state off half-filling (which is applicable to Yb^Asz, with rif ~ 0.1). The actual 
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lattice coherence scale, T co h, can be vastly different from the single-impurity Kondo scale, 
T K [14]; The c-fermion selfenergy can be computed from the d = oo approach [15]; in the 
low-T regime, it takes on the simple form, 

ImE c (u) = ^-[u 2 + 7t 2 T 2 ]x(l - x) (8) 

where D is the free bandwidth of the As (2p band(s)) and x = rif. This immediately tells 
us that the resistivity follows the T 2 law characteristic of FL behavior. The coefficient of 
the T 2 term is given by A = and is large for U » D. Notice that external pressure 
would increase D, reducing the coefficient A, as indeed observed. 

One might wonder why the shorter chain that gives rise to the observed ID magnetism 
does not seem to affect the transport properties. To see how this can come about, recall 
that in the spin-charge separation scenario, the holon and spinon velocities are Vh — xtf and 
v s ~ I e ff. These are much smaller than the Fermi velocity of the 2p carriers, as deduced 
from the bandstructure calculation [6]. Given that vp » Vh,v s , the scattering process 
involving the 2p hole and a spinon or a holon is unfavorable from the viewpoint of energetics 
(this process would dominate if v s < vf < Vh). Thus the heavy fermion FL characteristics 
involves, in our picture, the scattering of the light 2p holes off the (almost localized on 
the time scale of the p fermion hopping) localized / spins on the other three chains in 
Yb^As^. In this context, we mention that it has been proposed in [17] that scattering of 
the As p holes off the spinon excitations of the squeezed spin chain gives a resistivity going 
like T 2 . If this were indeed true, one would expect the T-dependence of the resistivity to 
change upon the application of an external magnetic field (which induces a spin gap) in 
contradiction with experimental evidence. The small magnitude of the magnetoresistance 
and the p(T, H) ~ A(H)T 2 behavior seen experimentally is thus in contradiction to the 
proposal in [17], but is completely consistent with the mechanism proposed here. In any 
case, it is highly improbable that scattering of light fermions off the spinons (holons) could 
yield such a large T 2 term in the resistivity. 

If the picture proposed above is correct, a more perfect charge ordering of the Yb f holes 
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along chains should disfavor the onset of heavy fermion FL formation. This connection is 
indeed borne out by experiments on F64P3 where a more perfect CO state is accompanied 
by insulating behavior with a larger 7 co-efficient. The more perfect CO state is a signature 
of decreased effective / — d — p hybridization, a corresponding decrease in the tendency 
to itinerance, and is consistent with the insulating resistivity observed for this material. 
A more perfect CO state also increases the number of spins on the shorter chain, and the 
higher 7-coefficient is easily understood as arising from the spinon excitations of the squeezed 
Heisenberg chain, as shown before. 

In conclusion, we have proposed that the apparently strange features in the low-T metallic 
state in Yb^Asz, i.e, a heavy fermion metal coexisting with ID quantum antiferromagnetism, 
are readily understood to arise from a single observation: the charge-order arising at higher 
T is not complete, and a small fraction of the S — 1/2 moments on Yb reside on the 
three longer chains. The ID antiferromagnetism is naturally understood in terms of the 
squeezed Heisenberg chain, while the heavy fermion FL properties at low-T result from the 
compensation of the Yb moments on the other three chains (the problem can be treated as 
three-dimensional) by the low density of carriers (coming from the As 2p bands). 
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FIG. 1. Magnetic susceptibility \ a t low temperature T for the spin S = 1/2 antiferromagnetic 
uniform Heisenberg chain of length N. 
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